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Phase transitions of first and second order can easily be distinguished in small systems in the 
microcanonical ensemble. Configurations of phase coexistence, which are suppressed in the conven- 
tional canonical formulation, carry important information about the main characteristics of first 
. order phase transitions like the transition temperature, the latent heat, and the interphase surface 

ON ' tension. The characteristic backbending of the microcanonical caloric equation of state T(E) (not 

to be confused with the well known Van der Waals loops in ordinary thermodynamics or mean field 
approximations) leading to a negative specific heat is intimately linked to the interphase surface 
'-^ ' entropy. 
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Microcanonical thermodynamics describes the dependence of the volume fl^ of the N-body phase-space of an inter- 
acting many-body system on the globally conserved quantities like the total energy, momentum, angular momentum, 
mass, charge etc.. This is the fundamental starting point of any statistical definition of thermodynamics. By Laplace 
transform of Qn(E) from the extensive quantities like the energy to intensive ones like inverse temperature (f3 = 1/T) 
I one obtains the more familiar Gibb's canonical partition function Z(JS). For systems interacting by short range two- 
body forces (with hard cores) both formulations are identical in the thermodynamic limit of infinitely many particles 
at the same particle density |IJ]. However, the two formulations are different for finite systems and more essentially 
are different even in their physical content for systems with long range forces like gravity ^ or Coulomb dominated 
\ systems see for example ||. This will be discussed in a forthcoming paper, 
"p-; ■ One of the most interesting phenomena in thermodynamics are phase transitions. Computer simulations of simple 
models give a deep insight into the mechanism. Naturally these calculations can only be performed for small systems. 
The main characteristics like the transition temperature T tr , the specific latent heat qi at , and the interphase surface 
tension a sur f are extrapolated to the infinite system. In the thermodynamic limit it should not matter whether 
the calculations for the finite systems are performed canonically or microcanonically. However, we will show for 
the two-dimensional, q = 10 states, Potts model, which has a clear phase transition of first order from ordered to 
O ■ disordered spins, the extrapolations converge faster when started from the microcanonical-ensemble. Moreover, our 
investigation will clarify some of the most striking features of finite micro-ensembles, the occurence of the S-bend in the 
microcanonical caloric equation of state T(E) and the occurence of a negative specific heat at a phase transition of first 
order, see e.g.: refs. {Here we do not discuss the negative specific heat in thermodynamically unstable systems 

with long-range interaction as were discussed by Thirring The systems we discuss here are thermodynamically 
stable with short-range interactions. They behave properly in the canonical ensemble and have a positive specific 
heat there and in the thermodynamic limit.} 

The Hamiltonian of the Potts model is defined as 

H = E'{ 1 - , w- (!) 

i<j 

on a two dimensional lattice (here with periodic boundary conditions) of N = L * L spins with q — 10 possible values 
(components) . The sum is over pairs of nearest neighbour lattice points only and <7j is the spin state at the i-th lattice 
point. For a microcanonical ensemble we have to calculate the partition sum over all possible different configurations 
v with the same total energy E : 

E U ,E • (2) 

v 

This is done by using the " Ergodic Microcanonical Metropolis Monte Carlo Algorithm" (MMMC) HJnJ . [The main 
idea here is to sample the energy in a narrow band E — 4 < E{ < i? + 4 by Metropolis sampling. From the 
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slope of ln(P(Ei))\E, in this energy-band one can determine f3 m icro(E)]. fijv(-E), the specific entropy s(e), and the 
thermodynamic temperature are related by (we put Boltzmann's constant k = 1) : 

Ns(e) = ln{Sl N (E = N * e)} (3) 

a ( \ — (A\ 

Pmicro\€ ) — \*) 

T ^ = T^~7e)' (5) 

Hmicro V c J 

The connection to the canonical (Gibbs-) free energy F(f3) is via the Laplace transform in saddle-point approximation: 

Z((3) = / Q N (Ne) e~ pife Nde (6) 



o N[s(.e)-0e] T ^27tN/c(e) (7) 
^)k = /3 (8) 

F{fi) = -T*ln{Z{fi)} (9) 



c(e) 



de 2 (9«V 
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Plotting s(e) and f3 m icro(e) v.s. e we see a convex intruder in s(e) (fig.l), where s(e) is reduced compared to its 
concave hull, the straight tangent line {~ the canonical s(e)} by As sur f. For an infinite system van Hove's theorem 



■13 forbids the existence of any convex part in s(e), because then the system gains entropy if it would divide 
spontaneously into two equal pieces one with ordered spins and with entropy s\ = s(ei), the other with disordered spins 
and entropy S3 = s(e3). Both pieces together would have the larger entropy s = (si +ss)/2 > S2 — s(e% = (ei +e^)/2). 
This is very nicely discussed by Huller in several recent articles e.g. |Q. 

This argument does not apply to a finite system as the new surface of phase-separation reduces the specific entropy 
by Asph-sep which is proportional to the number of lattice spins being fixed in the dividing surface i.e. As p h- S ep °c 
L/N oc L^ 1 . With growing size of the system, L — > 00, f3 m i cro (e) approaches the horizontal line /3 C an(e) = Ptr, 
determined by the Maxwell construction of equal areas Ai = A r on the left side below (3 tr and on the right above it, 
lower half of Fig.l.. On the other hand the surface entropy As sur f is equal to one of these areas A, half of the shaded 
area under the S-like oscillation of l3 m i cro {e) between the curve and the line (3t r ■ (3 m i C ro{t) can directly be calculated 
by MMMC so also the area A = As sur f. As the Potts model has only nearest neighbour interactions the interphase 
surface tension cr sur f — T * As sur f consists of an entropy part only. As the entropy of the special configuration with 
planar phase-separation can only be < s(e) we have As sur f < As p h~ S ep = 0(1/ L). We make the conjecture that the 
area N * A under N * s(e) scales like the cross-section of the lattice (oc L). 

The intimate link of A — As sur f to the surface tension can also be seen from the way how the surface tension 
is calculated in conventional canonical Monte Carlo simulations e.g. |]l5|-]l7|: Here the probability -Fr(e) to find the 
system with temperature T at an energy e shows at the transition temperature two well separated maxima with a 
deep valley in between. The ratio of the minimum to the maximum is by m| equal to: 

P fr (min) 

N * As SU rf = —ln{ - , A (11) 

\l Ptr (maxVjPtr (max2) 

Fig. 2 shows L * As sur / vs. the inverse lattice size 1/L in comparison with the surface tension determined 
by the multi-canonical method [ fL8[ . (Actually the quantity shown is twice the surface tension as we use periodic 
boundary conditions and consequently two cuts have to be made to separate the two phases.) Both results scale 
to the analytically known asymptotic limit [[l6]jl8]. We make an interesting observation here : The microcanonical 
surface tension depends much less on the size of the lattice than the canonical one. That means one may get its 
asymptotic value already at smaller sizes of the system than canonically. 

Insight into the possible mechanism leading to this weaker finite size scaling is obtained from the projected corre- 
lation function g^(dx) [jl9|| . 

g W(dx = \i x -j x \) = jYl - "Wt^ ~jy)}- ( 12 ) 



1y ,Jy 
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In the canonical ensemble it approaches 

g {1 \dx) oc cosh{(dx - L/2)/£ (1) } (13) 
£ (1) =£AA+(2</L) 2 (14) 

with £ the correlation length of the infinite lattice fl9|| . Figure 3 shows g^\dx) at an excitation energy e — 0.7 in the 
middle of the coexistence region. In the microcanonical ensemble g^ 1 ' does not need to be an exponential-like function 
and in fact it drops down faster in the region of coexistence (0.33 < e < 1.05) . In any case, the region of coexistence 
of both phases, gets suppressed for large TV by the Laplace transform eq.(^|). It does not contribute because in this 
region there is no stationary point, no solution of the equation (||) and all information contained in events showing 
the coexistence of both phases is lost in the canonical ensemble. This will more thoroughly be investigated in a 
forthcoming publication. 

Most of these calculations were done on different workstations. For the largest lattice (L = 100) we performed 
3.2 * 10 6 sweeps. To destroy correlations due to the previous energy step the spins in an arbitrarily chosen rectangular 
part of the lattice were first uniformly lifted one value higher up (cyclic) and the energy was restored by switching 
random spins of the lattice. This procedure has the advantage that the perturbation per latice point is oc 1/L and 
all correlations inside the rectangle are conserved. In between these macroscopic moves the spins were individually 
updated in sequence one after the other. As can be seen from fig.l the accuracy for L = 100 is not optimal but 
sufficient for the demonstration of the method. Due to bad computational conditions here at the HMI it was not yet 
possible to do better. 

An important lesson can be learned here : Following M. Fisher [^o| a phase transition of first order is distinguished 
from a transition of second order by a non-vanishing interphase surface tension. That means the microcanonical 
caloric equation of state for a phase transition of first order must show an S-oscillation and consequently a negative 
specific heat c = de/dT. This solves the outstanding problem of the backbending and of the negative specific heat at 
a phase transition of first order in the microcanonical ensembles discussed before by many authors e.g. ref. It 
also appears at the condensation phase transition of selfgraviting astrophysical bodies which again must be treated 
microcanonical as the long range of gravity does not allow for the transition to the thermodynamic limit {A 
selfgraviting system at constant energy behaves microcanonically quite differently from the same system at constant 
temperature (canonical). Whereas in the second case it can totally implode at low temperatures and convert all 
surplus energy to the heatbath, in the first case it can only partly implode because all energy must be stored in the 
noncollapsed rest of the system.} 

Figure 4 shows the three crossing points ei, £2, £3 of the Maxwell line (3t r with (3 m i cro {e) as function of the inverse 
lattice length 1/L. The difference qi a t — £3 — £1 at L = 00 is the specific latent heat of the transition, which is 
analytically known from Baxter |2lJ. Our results up to the lattice length L = 100 agree well when extrapolated 
linearly in 1/L to L = 00. Only the middle solution of (3 m i C ro(^2) = Ptr fails. However, this point is numerically as 
well as also physically the instability point. For the latent heat it is irrelevant. Again the finite size scaling of the 
microcanonical ti turns out to be very weak. It is rather surprising to find this simple scaling at such small lattices. 
As the theory of finite-size scaling for the microcanonical ensemble docs not yet exist we have to take this promising 
result with care. However, as was discussed by Hiiller Q this may be at no surprise, as much of the size dependence 
of the canonical ensemble comes from the trivial factor N in the exponent of the Laplace transform (^|) which has 
nothing to do with the physics contained in the specific entropy s(e). In a future publication we will discuss possible 
reasons for this weak sensitivity of the microcanonical Potts model to the finite size of the lattice and we will compare 
the microcanonical correlation functions with the canonical ones. 

We believe our findings with the q = 10-states Potts model are characteristic for all microcanonical ensembles. In 
fact the backbending of the caloric equation of state T(E) was found in many other finite system e.g. fragmenting 
nuclei ||, fragmenting atomic clusters Q, "melting" phase transitions in van der Waals clusters |||| . More work, 
however, must be done to understand what the "surface" entropy, the area A under the S-shape of (3 m icro(^), means 
in these cases which cannot be extrapolated to infinite systems. The scaling property of the microcanonical ensemble 
is yet unknown and must be examined. 

We want to thank the Fachbereich Physik of the Freie Universitat Berlin and the Institut fur Hochenergiephysik, 
DESY, Zeuthen for the excellent computing conditions. Only with their generous support we were able to perform 
these numerical calculations. 
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FIG. 1. a) Specific entropy s(e) = J Q e /3 m i cr o(e)de vs. the specific energy e for the 2-dim. Potts model with q — 10 on a 
100 * 100 lattice. In order to visualize the anomaly of the entropy the linear function a + be (a — s(0.2119), b = 1.4185) was 
subtracted. Because we use periodic boundary conditions one needs two cuts to separate the phases and the convex intruder 
is twice the surface-entropy. 

b) Inverse temperature f3 m icro{t) = 1/T(e) as directly calculated by MMMC 
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'Specific Surface Entropy" As*L = AS/L 
of 2dixn 1 O-stat cs Potts model 
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FIG. 2. Surface entropy 2As sur f evaluated as half the shaded area in figure 1 as function of the inverse lattice length 1/L 
compared to the same quantity from the multi-canonical method of Billoire . 
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Figure 4 

FIG. 4. The three solutions of /9 m icro(ei,2,3) = Ptr as function of the inverse lattice length 1/L. 
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